AD-A058  542  CENTER  FOR  NAVAL  ANALYSES  ARLINGTON  VA  F/G  20/13 

DIFFUSION  THEORY  OF  REACTION  RATES.  II.  ORNSTEIN  - UHLENBECK  AP— ETC(U) 
FEB  78  M MANGEL 

UNCLASSIFIED  CNA-PP-230  ML 

w HI  IHi^HiHH 

■ I ■ ■ ■ ■jJtE.il  UdH  ■■■■  I-'  -—4  I II 


UULTISimillgLllli 


55  0002S0, 


DIFFUSIONJTiEORY  OE.REACTION 
PRNSTEIN^-  UHLENBE'CK  APPROX 


Professional 


DISTRIBUTION  STATEMENT  A 

Approved  loi  public  release; 
Distribution  Unlimited 


D D C 

SGaEIMG 


The  ideas  expressed  in  this  paper  are  those  of  the 
author.  The  paper  does  not  necessarily  represent 
the  views  of  the  Center  for  Naval  Analyses. 


CENTER  FOR  NAVAL  ANALYSES 

1401  Wilson  Boulevard 
Arlington,  Virginia  22209 


Marc  Mangel* 
February,  1978 


*Center  for  Naval  Analyses  of  the  University  of  Rochester, 

1401  Wilson  Boulevard,  Arlington,  Virginia  22209 

* 1 


ffBHHMHMUBHMiniMlAi 


ABSTRACT 


The  diffusion  theory  of  reac- 
tion rates  is  considered  in  the 


i 


Ornstein — Uhlenbeck  approximation. 
The  rate  constant  is  calculated  by 
using  a modified  Kramers  (eigen- 
value) method  and  by  using  the 
expected  time  formulation.  The 
transmission  coefficient  is  calcu- 
lated. It  is  shown  how  correlation 
functions  can  be  calculated.  In 
the  last  section,  a model  problem 
is  considered. 
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DIFFUSION  THEORY  OF  REACTION  RATES,  II 
ORNSTEIN— UHLENBECK  APPROXIMATION 


SECTION  1.  INTRODUCTION 


In  the  previous  paper  (1),  denoted  by  I,  we  derived  the  dif- 
fusion theory  of  reaction  rates  (equations  2.16-2.29  of  I).  In 
the  Einstein-Smoluchowski  (ES)  limit  of  large  viscosity,  we  were 
able  to  calculate  transmission  coefficients,  rate  constants  and 
lifetimes  of  activated  species.  The  calculations  were  relatively 
straightforward  because  the  ES  limit  eliminates  one  of  the  inde- 
pendent variables  (the  velocity)  so  that  partial  differential 
equations  simplify  to  ordinary  differential  equation  or  to  gradient 
systems.  In  this  paper,  we  do  not  use  the  ES  limit,  but  consider 
the  full  OU  equations  of  I.  Unlike  in  the  analysis  of  I,  we  gain 
little  by  considering  one  dimensional  problems  in  the  OU  approxi- 
mation. The  multidimensional  OU  equations  are  (see  I)  : 


= £6ij  + Fi(x,i^  + vi 

9v  9v3  9V1  9X1  9V1 


(1.1) 


ft  = J>!“  - V1  ifij  - -ij  ((F1  - V^p) 

9v  9v3  9x  9v  \ / 


(1.2) 


. 2 - 


— u (x,  v,00)  = £6 


i 9t 


i3  + F1  + v1  ^ - v1 

9v19vD  9V1  9X1  9V1 


(1.3) 
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In  (1.1-1. 3),  repeated  indices  are  summed  from  1 to  n.  Namely 
(xfv)  is  a point  in  the  2n  dimensional  phase  space.  The  force 
F(x)  is  derived  from  a potential  F1  (x)  = -3V/3X1;  V(x)  is 
assumed  to  be  known. 

In  the  above  equations 


F (x)  = -W(x) 


(1.3a) 


where  V(x)  is  the  potential  surface  sketched  in  Fig.  1.  The 
function  u(x,v,t)  is  defined  by 


u(x,v,t)  = E | f^x(t),v(t)^  | x ( 0 ) = x,  v(0)  = vj 


(1.4) 


where  f(z,z  ) is  an  integrable  function.  By  appropriate  choices 
of  f,  we  can  construct  the  moments  of  (x,v) . The  density  for 
(x(t),v(t))  is  p (x, v,t) : 

p (x,  v,  t)  dxdv  = Pr  |x(t)6(x,x  + dx)  , v(t)€(v,v  + dv)|  (1.5) 
Finally  t(x,v)  is  an  expected  time 


■v)  = I * it  <x'v- 


t)dt 


(1.6) 


In  this  paper,  will  will  use  equations  (1.1-1. 3)  to  calculate 
reaction  rates  and  transmission  coefficients.  Our  methods  use 
the  asymptotic  techniques  developed  in  (2)  and  (3).  Since  these 
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techniques  are  quite  involved,  we  will  treat  the  problems  in  an 
order  slightly  different  from  that  of  I. 

In  addition  to  more  complicated  techniques,  the  formulation 
of  the  problems  is  more  delicate.  To  illustrate  some  aspects  of  the 
delicacy,  we  will  consider  the  phase  plane  in  both  OU  and  ES 
cases.  (i.e.  the  special  case  x€R  ).  The  phase  plane  corresponds 
to  the  deterministic  equations  obtained  by  averaging  the  equations 
(I.  2.17,  18,20.)  These  give 


& = F(X) 


= F(x) 


(1.7) 


(1.8) 


The  ES  "phase  plane"  is  just  a line,  sketched  in  Fig.  2.  In  this 

A 

case,  formulating  problems  in  terms  of  rates  to  cross  x,  or  times 

/\ 

to  reach  x is  perfectly  sensible.  The  OU  phase  plane  is  partially 
sketched  in  Fig.  3 (in  §6  we  give  more  detail) . There  is  a 

A 

trajectory  S that  plays  the  role  of  x in  the  ES  case:  S 

separates  domains  of  attraction  of  (Xq,0)  and  (x^,0)  (double 
minimum  case)  or  (°°,0)  (Lennard-Jones  case)  (see  I).  This  tra- 
jectery  is  the  deterministic  separatrix.  In  the  following,  we 
shall  use  the  separatrix  quite  a bit. 


i 


FIG.  2:  THE  ES  PHASE  SPACE. 

A)  DOUBLE  MINIMUM  POTENTIAL  SURFACE. 

B)  LENARD  JONES  POTENTIAL  SURFACE. 
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For  co  nceptual  ease  we  connect  the  phase  plane  with  the 
phase  space,  by  assuming  that  the  potential  function  V(x)  is 
quasi-separable . For  example,  if 

V(x)  = V(x2,...xn)  - jfx1) 3 + ax1,  (1.9) 

then  the  (x1^1)  phase  plane  will  be  similar  to  the  one  sketched 
in  Fig.  3.  The  full  phase  space  will  be  the  direct  product  of  the 
(x^v1)  phase  plane  and  the  2n  - 2 dimensional  (x2, . . .xn,v2,  . . . 
vn)  phase  space.  Hence,  we  are  assuming  that  the  reaction  co- 
ordinate is  clearly  defined  (see,  e.g.  (4)). 

In  §2,  we  construct  the  reaction  rate  constant  by  calculating 
the  lowest  eigenvalue  of  the  forward  equation  (1.2).  We  assume 

that  the  potential  is  of  the  form  (1.9)  and  that  the  reaction  pro- 

2 n 

ceeds  along  a path  in  which  (x  ,...x  ) have  equilibrium  values. 
This  assumption  is  made  for  notational  convenience.  In  §3,  we 
construct  the  reaction  rate  constant,  under  the  same  assumption 
about  the  potential,  using  the  expected  time  formulation.  In  §4 
we  calculate  the  transmission  coefficient.  In  §5,  we  illustrate 
how  correlation  functions,  and  thus  spectra,  can  be  found  by  using 
the  "ray  method".  Finally,  in  §6,  we  consider  a model  problem 
that  illustrates  some  of  the  ideas  in  §2-5. 
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SECTION  2 

REACTION  RATE  CONSTANT  BY  MODIFIED  KRAMERS  METHOD 

In  this  section,  we  use  the  Miller  (5)  — Ludwig  (2)  theory 
to  compute  the  rate  at  which  particles  pass  from  the  well  at 

A 

(Fig.  1)  to  the  barrier  at  x.  We  seek  a solution  of  (1.2)  in  the 
form 


yp(x,v,t)  = 


£ o (x , v ) e Anfc 
n n 


(2.1) 


In  §2-6,  we  make  the  following  assumptions  about  the  reaction. 

They  are  made  to  ease  the  conceptualization  of  the  ideas  presented 

here  and  can  easily  be  relaxed.  First  we  assume  that  V(x)  has 

2 n 

the  form  indicated  in  (1.9).  Next,  we  assume  that  (x  ,...x  , 
v2 , . . . vn)  have  relaxed  to  their  equilibrium  values  in  a time  tre^ 
much  smaller  than  the  first  passage  time  t^  considered  below. 

Thus,  we  reduce  (1.1-1. 3)  to  equations  involving  only  x^  and  v^, 
which  we  denote  by  (x,v) . This  type  of  assumption  was  used  by 
Slater  (6)  and  justified  by  Rice  (7) . In  I we  showed  that  the 
lowest  eigenvalue,  X^,  is  asymptotic  to  the  rate  at  which  molecules 
pass  from  Xq  to  x^.  We  are  thus  lead  to  consider  the  problem 

(2.2) 


32o 


-Xo  = e Z-Z  - v|£  - | [ F (x)  - v)o 

3V*  x 


l 


o(x,v)  = 0 


for 

(x,  v)£S 

lim 

2 

v 

Q 

x" 

< 

II 

o 

We  now  consider  the  semi-infinite 

domain 

in 

(x,v) 

consisting 

of  all  points  to 

the  left  of 

S = 

3D. 

(2.2)  over 

D and  obtain 

-X  = 

- — ( 
3v  \ 

(f  (x)  - 

v]qN 
1 ) 

jdxdv 

| Jo  (x,v)dxdv 
D 


using  Gauss'  theorem,  we  obtain: 


J Jo  (x,v)dxdv 
D 


In  (2.6),  t is  the  tangential  coordinate  along 


v (t) 

satisfy 

Sf-’ 

dv 

d^  * 

F (x) 

with 

A 

x (°°)  = X, 

V (oo) 

= 0. 

(2.3) 


space  D 
We  integrate 


(2.4) 


(2.5) 

S and  x (t ) , 

(2.6) 
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Following  the  argument  in  I,  we  now  set 


o (x,  v) 


+ 


+ 0(e) 


(2.8) 


Then  using  Laplace's  method  on  the  integrals  in  the  numerator  and 
denominator  in  (2.6),  we  obtain 


X ~ 


2 

TT 


V 


V"  (x0)  | V"  (X)  | 


- V ( x ) / e V (x  ) / e 
e e 


(2.9) 


the  result  (I,  5.16).  It  is  interesting  that  the  results  (2.9)  and 
(I,  5.16)  agree  exactly.  They  agree  for  the  following  reason. 

In  the  use  of  the  above  asymptotic  method,  we  implicitly  assume  a 
rapid  relaxation  of  the  velocity,  so  that  the  particle  crosses  at 

A 

x with  v inf intesimally  small.  Consequently,  we  make  an 
implicit  ES  assumption,  leading  to  (2.9). 
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SECTION  3 

REACTION  RATE  CONSTANT  BY  EXPECTED  TIME 


The  analysis  in  the  previous  section  implicitly  assumes  that 


the  barrier  V(x)  - V(Xg)  is  large,  a point  first  noted  by 


Kramers.  (8)  For  many  reactions,  e.g.  catalyzations,  this  is  not 
the  case.  Hence,  the  procedure  used  in  section  2 breaks  down. 

An  alternative  calculation  of  the  rate  constant  uses  the 
expected  time  formulation  (see  I,  §5.2).  We  consider 


-1  = e ift  + F(x)  it  it  _ v at 

^7  ' ; av  ax  av 


(3.1) 


with  boundary  conditions 


t (x , v ) = 0 
(x,v)€S 


lim  t(x,v)  < oo 
2 


lim  t(x,v)  <® 

x->  -oo 


(3.2) 


As  defined  in  (3.2)  t(x,v)  is  the  expected  time  to  reach  S 


. 'X,  % 

conditioned  on  x(0'  = x,  v(0)  = v. 


Let  pQ(x,0)  be  the  initial  distribution  of  molecules;  it 


will  be  concentrated  at  (Xq,0) 


We  define  the  rate  constant  by 


k = l/<t (x, v) >, 


(3.3) 


11 


where 


<t(x,v)>  = JCp  (x,v)t(x,v)dxdv. 

I-N 


(3.4) 


In  (3.4),  D is  the  domain  consisting  of  all  points  to  the  left  of 
the  deterministic  separatrix  S. 

To  construct  t(x,v),  we  seek  a solution  of  the  form  (3) 


t (x , v ) = g (x , v,  e ) B (^(x,  v) /e1//3  , ct/e2/3 , l/e1/3 , y) 

+ e2/3h(x,v,£)B' (x,v)/e1/3/a/e2/3,l/e1/3. 


Y) 


+ k (x , v,  e) 


(3.5) 


In  (3.5)  i Mx,v),  g (x,  v)  , h(x,v),  k(x,v,e),  a and  y are  to  be 
determined.  We  assume  that  they  have  asymptotic  expansions  (3) 

g (x,  v,  e)  = X)gn(x,v)en  h(x,v,e)  = £hn(x,v)en 


k (x,  v,  e)  = 2Zkn(x,v)en  a = X/i  e11 

n 


(3.6) 


En 

^ne 


The  function  B(z,8,l,y)  satifies 


d2B  .2  dB  . 

— 1 = " (z  “ 6)  d7  “ 1 + 

dz  az 


(3.7) 
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We  shall,  in  fact,  only  consider  the  first  term  here.  Higher  order 
terms  are  treated  in  an  analogous  fashion.  When  derivatives  are 


evaluated,  (3.7)  is  used  to  replace  B"  by 


\ B'  - 1/e1/3  • 


+ yi^/e1/3.  After  derivatives  are  evaluated  and  substituted  into 
(3.1),  terms  are  collected  according  to  powers  of  e.  The  leading 
term  vanishes  if  the  following  equations  are  satisfied: 


(3.9) 


0 ( e-1/3B ' ) : ^F(x)  - v^v  + V ^ ^ - a^)  = 0 (3 

0(e°B):  ^F(x)  - v^g®  + vg®  = 0 (3 

0 ( e° ) : ^F(x)  - v^k°  + k°  + ^g°(-l  + yQip)  = -1  (3.10) 

0 ( e2/3B ' ) : ^F(x)  - v^h°  + vh°  - h°^v(4-2  - ctQ))v  (3. 

- (*P2  ~ aQ)h°<pv  - a^^g0  - (<P2  - ao)h0)  + g°Kv  = 0 

We  first  consider  (3.8).  In  I,  we  showed  that  the  classical 
(deterministic)  equations  for  x,  v are 


(3.8) 


0 (g2/3B 


(3.11) 


= F (x)  - v 


(3.12) 


These  correspond  to  a set  of  "deterministic"  trajectories  in  the 
(x,v)  phase  plane;  parametrized  by  x(0)  = x and  v(0)  = v.  The 

A 

points  (Xq,0)  and  (x,0)  are,  respectively,  a metastable  and 
unstable  equilibrium  of  the  system.  The  system  is  not  a Hamiltonian 
system,  due  to  the  -v  term  in  the  second  equation  in  (3.12). 

Also,  we  are  interested  in  an  explicitly  non-equilibrium  phenomenon. 
Consequently,  we  require  ^ l*  0 at  either  of  the  rest 

points.  At  these  points,  (3.8)  becomes 


4>v('P2  ~ aQ)  = 0 (3.13) 

We  set  (3) 

<MX0 , 0)  = ~ <Mx,0)  = y[a^  (3.14) 

a 

The  derivatives  ^x'^v  at  (xQ,0)  and  can  be  determined  by 

differentiating  (3.8)  with  respect  to  x or  v and  evaluating 
the  result  at  the  rest  points.  Once  these  are  known,  the  parameter 
otg  can  be  calculated  by  the  method  of  characteristics.  Once  otg  is  known 
ip(x,v ) can  be  completely  determined  by  the  method  of  characteristics  (3,9). 
An  argument  using  Hamilton- Jacobi  theory  (3)  shows  that 


<()  = on  s>  Hence  the  separatrix  is  a contour  of  ip.  If 

we  set 


(3.15) 
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then  the  first  two  terms  on  (3.5)  vanish  identically  on  S. 

Equation  (3.9)  incidates  that  g°  is  constant  on  trajectories. 

y\ 

At  (Xq,0)  and  (x,0),  equation  (3.10)  becomes 


*v2(xo'0)g0(-1  ■ Yo \Po  } = _1 


(x'°)9° (~i  + y0  >/%  ) = -1 


(3.16) 


(3.17) 


These  are  two  equations  for  the  constants  g and  Next,  we 

set  )Cq(x,v)  5 0 on  S and  solve  (3.10)  by  the  method  of  charac- 
teristics (9) . With  the  above  choice  of  initial  data,  t(x,v) 
vanishes  identically  on  s. 

On  s,  equation  (3.11)  becomes 

0 


dh  0 j 2 0.  .2  0 

dt"  “2h  V a0  K = 9 Kv  + aiK9 


(3.18) 


At  (x,0)  we  obtain 


0 - ■(g0^vv  + “l^2^0) 

hu  (x, 0)  = — 


(3.19) 


Then  h can  be  determined  on  S,  and  in  the  entire  phase  plane 
by  the  method  of  characteristics.  The  parameter  is  determined 

in  a manner  analogous  to  the  determination  of  aQ. 


15 
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We  have  thus  given  a complete  characterization  of  the  first 


term: 


t (x, v)  = g°B  + e2/3h°B'  + k°  + 0(e) 


(3.20) 


With  this  construction,  we  can  obtain  the  rate  at  which  particles 
reach  S from  (x,v),  as  given  by  (3.3). 
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SECTION  4.  TRANSMISSION  COEFFICIENT 


In  this  section,  we  give  an  asymptotic  technique  for  the 

calculation  of  the  transmission  coefficient.  As  before,  rather 

/\ 

than  using  (x,0)  as  a reference  point,  we  use  the  separatrix  S 
(Fig.  4).  We  surround  S by  a tube  (S^,S2).  Let 


u(x)  = Pr 


= x,  v(0) 


|^x  (t)  ,v(t)^ 
v(0)  = vj 


crosses  S2  before  S^jx(O) 


(4.1) 


If  si»S2  are  sufficiently  far  apart,  then  it  is  clear  that  u(x,v) 
represents  the  probability  that  a molecule  with  initial  reaction 
coordinate  x and  reaction  coordinate  velocity  v becomes  a 
product  rather  than  returning  to  a reactant.  The  function  u(x,v) 
will  satisfy 


» - ^ + * ('<«>  - ’)£ 


(4.2) 


with  boundary  conditions 


u(x,v)  = 0 (x,v)6  S. 


(4.3) 


u (x , v)  = 1 (x,v)€  S. 


(4.4) 


I 


A formal  asymptotic  solution  of  (4.2)  is  (3) 

u(x,v)  = g(x,v,e)E  (^>  (x  ,v)  / + e1//2  h (x , v,  e ) E (x,  v)  / yJT  ^ 

(4.5) 

with  ip  (x,v)  ,g(x,v,e)  and  h(x,v.e)  to  be  determined.  We  assume 
that  g and  h have  asymptotic  expansions 


g (x,  v,  e)  = y^engn  (x,v) 
h (x,v, e)  = y^engn (x.v) 

In  (4.5),  E(z)  is  the  error  function: 


(4.6) 

(4.7) 


d2E  _ dE 

'Z  K 


z0  - z - Z1 


(4.8) 


Following  a procedure  analogous  to  the  one  in  the  preceding  section, 
we  obtain  equations  for  ip,g®  and  h^.  (Higher  order  terms  are  treated 
treated  in  an  analogous  fashion)..  These  equations  are: 


v*x  + 1 

^F(x) 

V 

- v^v  pip2  = 0 

(4.9) 

vg°x  + 

^F(x) 

i 

Jfl 
< o 

II 

o 

(4.10) 

vh°  + 

X 

^F(x) 

- v | h®  + g%  - h%  ^ 

/ v vyvy 

(4.11) 

0 ,0 
+ g ip  v - h 
v 


((Wv>v)  = 
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A 

We  require  tpx  ^ 0 , ? 0 at  (x,0)  . Then  (4.9)  indicates 

/v 

that  !p(x,0)  = 0.  An  argument  using  Hamilton- Jacobi  theory  shows 
that  ip  H 0 on  S.  When  (4.9)  is  differentiated  with  respect  to 
x,v  and  evaluated  on  S,  we  obtain 


jT  ip  + F (x)tl/  - tp  ip  = 0 
dt  rx  x ' 


xTv 


(4.12) 


^ + ip  -tp  - ip3  = o 

dt  v yx  yv 


(4.13) 


In  obtaining  (4.12,  13)  we  used 


at  \ = v\x  + 


^F(x)  - v)t 


(4.14) 


In  light  of  (4.14),  at  (x,0)  (4.12,  13)  become 


Fx(x)lpv  “ = 0 


(4.15) 


ip  - ip  - ip"  = 0 
rx  yv  Yv 


(4.16) 


We  denote  the  solution  of  these  equations  by  tp  ,tpv.  Then 

X * 

(4.12,  13)  are  solved,  subject  to 


tPxM  = *x 


ip  (oo)  = ip 

YV  yv 


(4.17) 


Thus,  tp  and  tp  are  known  along  S.  Consequently,  i|j(x,v)  is 
known  in  the  entire  plane  since  (4.9)  can  be  solved  by  the  method 
of  characteristics,  with  initial  data  on  S (3,9). 
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Let  denote  the  minimum  value  of  ip  on  S ^ and  ^ 

the  maximum  value  of  ip  on  S . Equation  (4.10)  indicates  that 
g°  is  constant  on  trajectories.  Following  (3),  we  set  zQ  in  (4.8) 
equal  to  ipj./  and 


90  = E«’ll/ 


V^>] 


(4.18) 


If  S and  S2  happen  to  be  level  curves  of  ip,  then  the  choice 
(4.18)  insures  u(x,v)  ho  on  and  u(x,v)  = l on 

Otherwise,  u is  exponentially  small  on  and  1 - u is  ex- 

ponentially small  on  S2. 

Finally,  we  consider  the  equation  (4.11)  for  h°.  On  S, 
it  becomes 


dh  .0,2  0, 

-tt — - h ip  = -g  \p 
dt  rv  ^ vv 


(4.19) 


At  (x,0)  equation  (4.19)  gives 


,0  0-  /T  _ .0,  . 
h = g ^vv/^v  = h (°°) 


(4.20) 


where  ip  = t|/vv(x,0).  Then,  on  S,h  is  given  by 


h ( t)  = 


r°°  o r rs  2 1 o 

9 ^vv(s)  exp  -J  v|>v(t’ ) dt’  ds  + h (°°) 
exp  - J ^v(t* ) 2dt • 


(4.21) 


i 
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SECTION  5 

CORRELATION  FUNCTIONS  AND  SPECTRA 

We  now  show  how  to  calculate  correlation  functions  and 
their  Fourier  transforms.  We  define 

4>  xx  (t)  = E jx(t  + x)x(t)| 

= E jx(t  + T)v(t)j=  E jx(t)v(t  + t)| 

{’ 

The  spectra  are  obtained  by  Fourier  transforms,  e.g. 

Ikl  (a>)  = J (cos  ^t)  ^ (T)dx 


(t)  = E { v ( t +x)v  (t)| 


(5.1) 

(5.2) 

(5.3) 


(5.4) 


Each  of  the  spectra  is  calculated  in  a similar  fashion,  so  that  we 
only  consider  <|>xx(x),  which  for  simplicity  we  denote  by  4>(t). 


If 

4> 


xv 


(t)  = E |x(t  +T ) x ( t ) |x(t)  = x,v(t)  = v| 


(5.5) 


is  the  conditional  correlation  function  and  p(x,v,t)  is  the  density 


,0. 


for  (x,v)  at  time  t,  then 

<f>  (t) 


= JJ  <J>XV  (x)  p (x,v  ,t)dxdv 


(5.6) 


The  ray  method  can  be  used  to  satisfactorily  calculate  ^xv^T^' 
by  integrating  along  rays  that  emanate  from  the  point  (x,v) , (see 
I,  §6  and  Fig.  5) . 
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SECTION  6.  AN  EXAMPLE 


In  this  section  we  consider  the  dynamical  equations 

dx  = vdt  (6. 1 

2 

dv  = (x  - a - yv)dt  (6.2 

3 

The  potential  function  V(x,a>  = -l/3x  + ax  is  sketched  in  Fig. 

In  Figs.  7a-7e  we  sketch  the  (x,v)  phase  plane  as  a varies,  f 

constant  y.  As  a decreases  to  zero,  the  height  of  the  potential 

3 /2 

barrier  (Q  = 4/3a  ) also  decreases.  In  Fig.  8 we  sketch  con- 

tours of  the  transmission  coefficient  for  e = .01. 


] 


FIG.  6:  POTENTIAL  FUNCTION  FOR  THE  SYSTEM  CONSIDERED  IN 
SECTION  6. 
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